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Domain-wall fermions and chiral symmetries 

Sinya Aoki*^ 

^ Institute of Physics, University of Tsukuba, Tsukuba, Ibaraki 305-8571, Japan 

We investigate chiral properties of the domain-wall fermion (DWF) system. After a brief introduction for the 
DWF, we summarize the recent numerical results on the chiral properties of the domain-wall QCD (DWQCD), 
which seem mutually inconsistent. We next derive a formula which connects a chiral symmetry breaking term in 
the five dimensional DWF Ward-Takahashi identity with the four-dimensional hermitian Wilson-Dirac operator. 
Based on this formula, we propose a solution, which resolves the inconsistency among recent numerical data, and 
give a consistent picture of the chiral properties of the DWQCD. 


1. Introduction 

A suitable definition of the chiral symmetry has 
been a long standing problem in lattice field the¬ 
ories. Recently an ultimate solution to this prob¬ 
lem seems to appear in the form of the Ginsberg- 
Wilson relation [||j|]. Two explicit examples of 
the lattice fermion operators which satisfies the 
Ginsberg-Wilson relation have been found so far: 
One is the perfect lattice Dirac operator con¬ 
structed via the renormalization group transfor- 
mationS and the other is the overlap Dirac op¬ 
erator derived from the overlap formalism or 
from the domain-wall fermion(DWF)|^J^ in the 
limit of the infinite length of the 5th dimension. 
Since the explicit form is simpler for the latter, a 
lot of numerical investigations ||^-|Ti|] as well as an¬ 
alytic considerations have been carried out 

for the domain-wall fermion or the overlap Dirac 
fermion. 

In this report we review the recent numeri¬ 
cal investigations on the chiral properties of the 
DWF, which, however, bring puzzling results. To 
resolve these puzzles, we derive a formula which 
connects a chiral symmetry breaking term with 
the 4-dimensional hermitian Wilson-Dirac opera¬ 
tor. Using this formula and the eigenvalue distri¬ 
bution of the Wilson-Dirac operator, we propose 
a consistent interpretation of puzzling results. 

2. Domain-wall fermions 


tion|^ ?], which is equivalent to the 5-dimensional 
massive Wilson fermion with the free boundary 
condition in the 5-th direction, interacting with 
the 4-dimensional gauge field. Symbolically the 
action is given by 

Sf = +j^d5]ip 

+ + d'^]ilj + Mip-ip. (1) 


Note that the sign of the mass term M is opposite 
to the usual Wilson fermion. Since gauge fields 
are 4-dimensional, the 5-th coordinate s can be 
interpreted as “flavors”, so that the action can 
be rewritten in the 4-dimensional form: 


Sf — Un^fi'4’n+ii,s 


+ [mpr + m^Pl] y™,*. 


( 2 ) 


where the last term is considered as the flavor 
mixing “mass” term. 

Now we consider the fermion spectrum derived 
from this action for the free theory. In this case 
the action in the momentum space becomes 


Sf = J d^p[ipsn^ 

+ ^s{MPr + M^Pl)^^H 

where 


/ -W 1 
-W 


M = 


-W 


\ 


1 

-W j 


( 3 ) 


Throughout this report we employ a Shamir’s 
variant of the original domain-wall fermion ac- 
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Table 1 

Relation between M and number of poles 


16 X 32 XN 3 , RG, P=2.6, M=1.8 


condition 

Pa* 

# of poles 

0.03 

0 < M < 2 

( 0 , 0 , 0 , 0 ) 

1 


2 < M < 4 

(tt/o, 0 , 0 , 0 ) 

4 

0.02 

4 < M < 6 

(tt/o, 7r/a, 0,0) 

6 


6 < M < 8 

(tt/o, 7r/a, 7r/a, 0) 

4 

0.01 

8 < M < 10 

(tt/o, tt/o, tt/o, tt/o) 

1 



= 


( -W 

1 -w 
1 



-w 

V i -w J 

with W = 1 — M + ~ cos(p^a)]. 

By setting i’{p)s = Prur{p)s + PlUl{p)s, 
one observe that the solution to the equation 
MstUR{p)t = 0 {MltUL{p)t = 0) gives a right- 
handed (left-handed) zero mode. The solutions 
to these equations are given by 

URip)s = W{py~^CR ( 4 ) 

ULip)s = W{p)^‘'~^CL, (5) 

where Ns is the number of sites in the 5-th dimen¬ 
sion. In order to satisfy the boundary conditions 

URiNs) = W{p)^‘-^CR = 0 (6) 

ul{ 1) = W{p)^^-^CL = 0 (7) 

for the above solutions, the Ns —>■ 00 limit is nec¬ 
essary, together with the condition that | lR(p) | < 
1 , which is equivalent to 

0 < M -f ^(cos(p^a) — 1 ) < 2 . ( 8 ) 

Since we have the kinetic term 7 ^ sin(p^a) for the 
fermion, the condition that ~ 0 or tt is re¬ 
quired for the existence of the massless pole in 
the continuum limit. Eq.( H) gives a condition for 
M and which is summarized in table |^. It is 
clear that the value of M (Wilson fermion mass) 
controls the number of massless fermion poles. 
For 0 < M < 2 we have one physical massless 
Dirac fermion. 

The form of zero-modes obtained above is too 
complicated in the coordinate space even for the 


Figure 1. as a function of m/ on a 16^ x 
32 X Ns {Ns = 4,10 and 20) at M = 1.8 and /3 = 
2.6 for the renormalization group (RG) improved 
gauge action. 

free theory to define the quark field as PrUr + 
PrUl- Instead we define the quark field as 

qn = Pr^u,! + PLi’n,Ns (9) 

Qn = ^upPl + i’n,NsPR- ( 10 ) 

This definition is reasonable since the right- 
handed zero mode has a peak at s = 1 and the 
left-handed one at s = Ns- 

According to this definition, the quark mass 
term is given by —mfqnqn, which leads to the 
non-zero pole mass of the fermion for 0 < M < 2 : 

mpoie = M(2-M)[m/ + (1-M)'^^]. (11) 

The formula tells us that rUpoie ~ (1 — y 0 
even at rrif = 0 for Ns y oo, corresponding to the 
fact that no solution to zero modes exist for finite 
Ns- On the other hand the fact that rupoie oc 
TO/ in the Ns —>■ oo limit strongly suggests that 
the chiral symmetry is realized for the DWF. The 
lattice QCD with the DWF is called the domain- 
wall QCD (DWQCD). 

As suggested by the behavior of the pole mass, 
the DWQCD satisfies the almost exact axial 
Ward-Takahashi identity, which is given by 

{{A^A\{t) + 2mfqj,T\{t) - 2jl{t,Ns/2)} O) 

+ {8\0) = 0 ( 12 ) 

where A 4 is an axial-vector current, O is an ar¬ 
bitrary operators which contain q and q only as 
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Figure 2. m^s as a function of Ng for the plaque¬ 
tte gauge action at a ~ 0.2 fm(circles) and 0.1 
fm(squares). 


the fermion variable, and j^{t, Ns/2), defined in 
the middle of the 5-th dimension, represents the 
explicit breaking term of the chiral symmetry. 

In perturbation theory it is proven that 
{j\{t,Ns/2)0) ^ 0 as Ns ^ oo for the non¬ 
singlet sector (6 7 ^ 0), while Ns/2)0) oc 

7 ^ 0 for the singlet sector (6 = 0). The 
latter one corresponds to the U^(l) anomaly, 
well-known in the continuum QCD. 

Again the presence of the exact axial Ward- 
Takahashi in the DWQCD suggests the existence 
of the chiral symmetry in the Ng —> 00 limit even 
for interacting theories. 


3. Numerical investigation of chiral prop¬ 
erties of the DWQCD 


In order to investigate the chiral properties 
of the DWQCD, we define the anomalous quark 
mass in axial Ward-Takahashi identity as follows. 


TO55 


(jg(t,jV,/2)-g75r°g(0)) 
t^oo (g75T“q(t) • q^5T°-q{0)) 


(13) 


Figure 3. Same as Fig.|| for the RG improved 
gauge action. 

renormalization group (RG) improved gauge ac¬ 
tion |pd| . 

The figure shows that mf dependence of is 
much milder than that of and is almost con¬ 
stant, so that the extrapolation of to m/ = 0 
becomes much easier. In addition 771^5 is very 
precise and sensitive to Ng , so that the Ng depen¬ 
dence can be extracted. Hereafter we exclusively 
employ m^s to investigate the chiral properties 
of the DWQCD, because of these advantages of 

TUqS. 

We have plotted niqs as a function of Ng, cal¬ 
culated on the CP-PACS for the plaquette gauge 
action in Fig.| and the RG improved gauge ac- 
tion(lower) in FigJ|, at coarse(a ~ 0.2 fm) and 
fine(a ~ 0.1 fm) lattice spacings[|l). 

At both lattice spacings is smaller for the 
RG improved action than for the plaquette ac¬ 
tion. At the coarse lattice spacing(a ~ 0.2 fm) 
it seems that C 7 ^ 0 for both actions. On the 
other hand, C = 0 for the RG action while C 7 ^ 0 
for the plaquette action at the fine lattice spacing 
(a ~ 0.1 fm). 


We expect that behaves as 

TO55 = C + aexp[-^Ng]+"fmf+ 0{m‘f).{l'i) 

The existence of the massless fermion requires 
C = 0. 

In Fig. we have plotted rriq^ as a function 
of mf with M = 1.8 on a x T x Ng = 
16^ X 32 X (4,10,20) lattice at P = 2.6 of the 


4. Phase structure of the Wilson fermion 
and the DWQCD 

In this section we try to understand the de¬ 
pendences of iriqs on the lattice spacings and the 
gauge actions. The existence of massless fermion 
of the DWQGD can be proven in the absence of 
zero eigenvalues of the hermitian Wilson-Dirac 
operator Hw = 'y^Dw{M), where Dw{FI) is the 
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Figure 4. The expected phase structure of the 
lattice QCD with the Wilson fermion in the (3 = 
-M plane, where tt = qi'y^r^q. 


PL,p=5.65,m„=1.7,12" PL,M-0.mo=1 -8,16h2 
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4-dimensional Wilson fermion operator with the 
fermion mass M and Dw{My = is 

satisfied. 

In the case of the 4-dimensional Wilson 
fermion, zero eigenvalues of JbDw leads to the 
spontaneous breaking of parity-flavor symme¬ 
tries, {qi'ydT^q) ^ 0, as shown belowQ. 


Figure 5. Monte Carlo time histories for the 
lowest six eigenvalues of obtained with the 
plaquette gauge action. 

gauge coupling constant [Q. The small value of (3 
corresponds to the coarse lattice spacing while the 
large value to the fine lattice. At = 0 (/3 = 00 ) 
no zero eigenvalues exist for 0 < M <2, which 
is the region where the domain-wall fermion has 
the zero modes. If we increase the gauge coupling, 
the allowed region of M for the massless DWF be¬ 
comes narrow and is shifted to larger values of M. 
Finally the allowed region disappears at /3 = /3c. 
This implies that massless fermions disappear at 
/3 < /3c for the (quenched) DWQCD. 

According to this interpretation it seems that 
/3c > 5.65 (a = 0.2 fm), 6.0(a = 0.1 fm) for the 
plaquette gauge action, while the numerical data 
indicate that 2.6 (a = 0.1 fm) > /3c > 2.2 (a = 0.2 
f^y for the RG improved gauge action. 


{qi^ 5 T^q) = — lim Tr 
— lim tr 


*75 T 


H^+O Dy\r + 

*75 *75 


ff-*+o [D\Y + i^ 5 H Dw — i'JbH 

■ r 1 1 

—i hm tr -^-^— 

H—>+o Hw + 'i'H Hw ~ iH 

—i lim > (A. 


( 1 

1 ^ 

V 

Xu-iHj 




1 


1 


X + iH X-iH 


= J dX pHw{X){-2'n:i)5{X) 

= -2tt pHw{^) 

where pHw ("^) is the density of the eigenvalues of 
Hw, which is defined by 

P**w(A) = ^<5(A„-A). 

n 

In Fig.^, we have drawn the expected phase 
structure of the lattice QCD with the Wilson 
fermion in the /3 = Q/g^-M plane, where is the 


5. Eigenvalues of Hw and PHwi^) 

(l&i this section we directly investigate the dis¬ 
tribution of small eigenvalues of Hw 

In Figs. H and ^ we plot Monte Carlo time his¬ 
tories for the lowest six eigenvalues of H^ for the 
plaquette and RG-improved actions |@. The fig¬ 
ures at a ~ 0.1 fm are shown at the same lattice 
size with the previous work of m^q. In each figure 
the left panel shows results for a ~ 0.2 fm and the 
right panel for a ~ 0.1 fm. Open squares show 
the minimum eigenvalue A^;^ and filled diamonds 
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Figure 6. The same as Fig. || obtained with the 
RG-improved gauge action. 


are the higher five eigenvalues. 

There is a clear trend that the minimum eigen¬ 
values become larger for smaller lattice spac- 
ings. Another interesting point is that the RG- 
improved action gives larger values than the 
plaquette action, which indicates that the RG- 
improved action has a better chiral behavior. 
The spectral density of Hw is defined by 

Pgw(A) = Jim ^ VJ(A^-A), (16) 

V—*oo o • 4 • K 


where the summation is over the eigenvalues of 
Hw- The density of zero-eigenvalues pHwi^)^ re¬ 
lated to the order parameter of the parity-flavor 
breaking, has been calculated by the accumula¬ 
tion method proposed in which is based on 
the relation 


1 

3-4-V 


E 1 = 

|A'|<A 



dAV(A') 


2p(0)A + O(Aj. 


(17) 


We note that, for the small-A expansion in (p^, 
analyticity of p(.^at the origin is assumed. 

The result of , where pHw (0) several /3’s 
have been plotted as a function of M, is given in 
Fig. g for the plaquette gauge action. Although 
PHwi^) becomes smaller for smaller a (; larger (3), 
no region without zero eigenvalues seems to ex¬ 
ist. This impression is indeed confirmed in Fig. g. 



1.0 1,5 2.0 


m 

Figure 7. pHwi^) a function of M at several 
values of (3 for the plaquette actionjl^. 


where phw (0) has been plotted as a function of 
a-i/2 at M = 1.2@. 

It has been speculated in [?] that the a depen¬ 
dence of phw (0) is well reproduced by the from 
that 

Pffw(0)-^exp[—(18) 

If this form is correct for all M, the region with 
Pffw(O) = 0 exist as long as a ^ 0. 

The similar result has been obtained for the RG 
improved action. 

These results immediately lead to the conclu¬ 
sion that there is no gap of the parity-flavor 
breaking phase in the g^-M plane for the Wilson 
fermion, so that the domain-wall QGD fails to re¬ 
alize the chiral symmetry at all a ^ 0. However 
this conclusion seems to contradict the numerical 
data for the RG improved gauge action at o = 0.1 
fm. Puzzles still remain. 


6. Theoretical understanding for the be¬ 
havior of mg 5 


In this section we try to understand the patho¬ 
logical behavior of as a function of Hg, in 
terms of the eigenvalues of Hw ■ 

We first derive the analytic expression for 
as follows. 


2mg5(J5 

(I -m/)2 


Ex.v 1 

!{H)-^H 

GW 

y)|2 


^GW 

y>|2 


(19) 
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10-5 -^^^^^^^^^^^^^^- 

1.5 2.0 2-5 3.0 

(aVa)”^ ^ 

Figure 8. PHwi^) ^ function of a“^/^ at M = 
1.2 for the plaquette actionfl^. 


where X = (x,a,a) represents site, spinor and 
color indices, 


m) 

T 




H, 


GW 


rp-Ns/2 _|_ rpNs/2 

cosh{HNs/2) = - - - 

Hw - 

1 

-=-75 

(1 -I- nif) -I- (1 - mf)j 5 tanh(FfiVs/2) 


For the derivation of the above formula, see . 

Based on the above formula can be ap¬ 
proximately estimated at to/ = 0 as follows. 


2'tTlq^Qj^ — 


N, 


all 


■ Nd 

Co ^ fiKY 

n—1 


(Nall - Nd) / dXp{X)fCX) 


-2 


1 


A = log 


Nall 

1-A 


1 + A 


coY.f{Xn)-^+Y.fCK)-Y20) 

c 

( 21 ) 


where Ac(A„) is the continuous (discrete) eigen¬ 
value of Nd is the number of discrete modes, 
Naii= L^xTxNaXA is the number of degrees of 
freedom, and cq is a weight factor for the discrete 
modes. See |17 for the detail of the estimation. 



Figure 9. The expected distribution of the eigen¬ 
values of Hw in the weak coupling region. 


The important point is that there are two types 
of eigenvalues, the continuous eigenvalues and the 
discrete eigenvalues. The eigenfunction for the 
former is the plain-wave like and spreads over 
whole space, while the one for the latter is expo¬ 
nentially localized. In the strong coupling region 
at P < Pc, the continuous eigenvalues can become 
zero, so that phw (0) Y 0- Therefore the DWQCD 
does not work as expected. (See [?] for more de¬ 
tails.) In the weak coupling region, on the other 
hand, the continuous eigenvalues has a minimum 
value (Ac)min; such that PHwY^c ^ (Ac)niin) — 0 
for the continuous eigenvalues. Instead of the 
continuous eigenvalues, discrete eigenvalues may 
become almost zero. The expected distribution 
of eigenvalues A in the weak coupling region is 
schematically drawn in Fig. 

It is likely that small (localized) discrete eigen¬ 
values are caused by dislocations. In particular, 
Q, the value of the topological charge in some 
configuration should move to the different value 
Q' during the simulation of the lattice QCD, so 
that the dislocation appears during the transi¬ 
tion from Q to Q'. If the dislocation appears, the 
Wilson-Dirac operator has an almost zero eigen¬ 
value at some value of M. Therefore it is very 
difficult to avoid small (localized) eigenvalues of 
i/iv and . 

From the approximated estimate for only 
near zero eigenvalues can contribute to it as Ns —>■ 
oo. Hence only small discrete eigenvalues are dan¬ 
gerous for the DWQCD in the weak coupling re¬ 
gion. The analysis in [|^ leads to the followings. 

In the Ng oo limit to^s vanish at finite vol- 
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ume (;Nau 00 ), since the probability of having 
exact zero eigenvalues is zero. 

If the infinite volume {]Naii —>■ 00 ) limit is taken 
before the Ng ^ 00 limit, there are two possibil¬ 
ities for the distribution of discrete eigenvalues. 
Unless the number of near zero eigenvalues in¬ 
creases linearly in Naii, they does not contribute 
to TO 55 , hence the DWQCD works fine. If it in¬ 
creases linearly in Naii, PHwi^) hence 
becomes non-zero, so that the DWQCD can not 
realize the chiral symmetry. 

The numerical data obtained so far seem to pre¬ 
fer the latter possibility. However more detailed 
investigations will be necessary for the definite 
conclusion. 

7. Conclusion 

The distribution of near zero eigenvalues of H\y 
plays a crucial role for the chiral properties of the 
DWQCD. We have pointed out that there are two 
types of eigenvalues, continuous and discrete. In 
the weak coupling region, only discrete eigenval¬ 
ues become relevant for mq^. However the effect 
of discrete eigenvalues to is rather small even 
if the number of near zero discrete eigenvalues in¬ 
creases linearly in the volume. We think that 
this has brought confusing numerical results for 
the Ng dependence of m^s. 

Our analysis suggests that the gap is closed in 
the phase structure of the Wilson fermion due to 
the near zero discrete eigenvalues, if the quenched 
approximation is employed. In the dynamical 
QCD simulation, the near zero discrete eigenval¬ 
ues are strongly suppressed by the fermion deter¬ 
minant, so that the gap opens and the expected 
phase structure is correct in the real lattice QCD. 
Since such a suppression of small discrete eigen¬ 
values is not known for the DWQCD, it is unclear 
how the dynamical DWQCD change the distribu¬ 
tion of the near zero discrete eigenvalues. 
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